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ABSTRACT 


A variable conductance heat pipe was designed and con- 
structed to study the performance characteristics of the 
gas loaded heat pipe with particular emphasis on the effect 
of gravity. Heat inputs were varied from 10 to 50 watts 
and the operating position of the heat pipe was varied from 
horizontal to vertical with the condenser above the evapora- 
tor.  Condenser temperature profiles are presented for 
(Pilous heat inputs and positions for the conventional and 
gas loaded heat pipe. The possibility of stratification in 
a horizontal heat pipe operating with a non-condensible 
gas that has a larger molecular weight than the working 
fluid is discussed based on observations presented. 
Padationally, a one dimensional analysrs™was ""Ormulaved 
and used to compare with the experimental results. The 
results show a good agreement between theory and experiment, 
ша сас те Спас the developed mathematical model adequately 


Ээ eribes the physical situation. 
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NOMENCLATURE 


A Cross sectional area of vapor space ft* 

А Cross sectional area of wick-wall ft* 

P сошот пас ой 

A Outside се астса о1 сопдјепсег ft^ 

SOC Gen 

D Diameter of vapor space (equals ci Е 

g Gravity ft/sec" 

h Heat transfer coefficient DU EX E 
ft 

Бос Latent heat of vaporization BTU/1bm 
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k Thermal conductivity ВОВЕ 

L Length of condenser section гү 

m Condensation rate lbm/ft^ sec 

Б Pressure lb/ft? 

Q Heat input | BTU/ft° sec 
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R. Inside radius of pipe IT 

T Local vapor space temperature ZR 

Щи Ambient temperature К 

V Voltage | volt 

W Mass ratio о” Ммогкіпр vapor 
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G Surface tension Dye 
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V Kinematic viscosity ft^/sec 





р Density 

u Dynamic viscosity 

015 Lennard-Jones force constant 

9, Collision integral based on Lennard- 
Jones potential function 

б Distance between wires in screen 
mesh wick 

SUBSCRIPTS 

a Working fluid 

b Non-condensible pas 

5 Stamlessrsteel 

L Liquid working fluid 
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ТЕ Dp СЛ ТИМ 


А. BACKGROUND 

Ehe heat pipe is a device capable of transferring large 
amounts of heat while operating at a nearly constant temper- 
ature. The major development of this device has occurred 
within the last decade although the original concept was 
proposed in 1944. This proposal, by R. S. Gaugler [Ref. 1] 
uu ne General Motors Corporation, was not acted upon at 
that time. A similar proposal by L. Trefethen [Ref. 2] was 
made to the General Electric Company in 1962 but was not 
investigated. The advent of the space program led to a 
renewed interest in the heat pipe in 1964 by G. M. Grover 
Seal |Ref. 3] of the Los Alamos Scientific Laboratory. 
The possibility of developing a reliable, efficient and 
Tiphtweignt heat transfer device for application in space 
equipment was too appealing to allow the concept to again 


be overlooked. 


ЕН HEAT PIPE OPERATION 

The sealed heat pipe is lined with a capillary wick 
EUCucture. The vessel is evacuated and partially filled 
mh a fluid referred to as the working fluid. During 
operation the working fluid is continuously evaporating 210) 
ме сестаол о: сре pipe and condensing throughout another 
section (see Figure 1). The condensate is returned to 


the evaporator section by means of capillary force acting 





threugh the wick structure. This change of phase process 
results in a nearly constant temperature throughout the 

heat pipe. The near isothermal operation of the heat pipe 
was reported in 1965 by K.F. Bainton [Ref. 4]. The iso- 
thermal characteristic of the heat pipe allows for a more 
efficient heat transfer device than that achieved with an 
equivalent solid metal rod. An effective thermal conduc- 
мању Гог а lithiumeheat pipe of 1,000”"to 10,000 times that 
of an equivalent metal rod was reported by workers at КСА 


ШІ 1967 [Ref. 5]. 





Heat pipe vapor flow diagram: 
figure 1 í 


The introduction of a non-condensable gas along with 
the working fluid alters the performance of the heat pipe 
appreciably. During operation the non-condensable gas is 


driven to the condenser along with the working vapor. 
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Unlike the vapor however the gas does not condense and 
thereby forms a high concentration of the non-condensable 
EE ina portion of the condenser. As the concentration of 
the non-condensable ass increases the amount of working 
fluid present and the condensation rate decrease. The 
mauserature of the pipe in this portien of Re condenser 
will therefore become less than in a portion vhere the поп- 
condensable gas is absent or is a very small amount of the 
concentration. Since the pressure of a saturated vapor is 
оше ссту proportional to its temperature, tHe amount of 
Working vapor present will decrease with the temperature. 
For a large enough concentration of non-condensable gas 

in a portion of the condenser, ambient temperature will 

be approached. The small amount of working fluid present 
forming the partial pressure of methanol will not condense, 
thereby reducing the heat transfer area and the conductance 
ОГ the heat pipe. 

If the evaporator temperature increases, the partial 
pressure of the working fluid increases thereby compressing 
tne non-condensable gas into a smaller volume which in 
turn increases the heat transfer area. The pipe's conduc- 
рапсе 15 thereby increased. This characteristic of the 
gas loaded heat pipe enables it to have the potential to 
perience smaller operating temperature variations than 
conventional heat pipe for similar changes in heat Input. 


Therefore, the process by which the non-condensable gas 


Ше) 





Tem 
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collects in the condenser has a direct bearing on the 
heat pipe's performance. 

Two analyses have been developed to give insight into 
pue behavior of the gas loaded heat pipe. The first, the 
flat front theory, was developed by W. Bienert [Ref. 6]. 
This theory assumes that during steady state operation of 
the gas loaded heat pipe the vapor and gas are separated 
cae condenser by a vertical That front. Тһе axial heat 
eonduction of the pipe is considered negligible in this 
Theory. The possibility of a pressure grag@itent being 
Би ес in the pipe 1s excluded and diffusion is not con- 
Sidered. With all of these assumptions an analysis of the 
heat pipe performance is straight forward and proves ade- 
quate in the evaluation of several heat pipe parameters. 
Specifically the volume of the gas reservo?r and the ratio 
of condenser to reservoir volume can be determined. Figure 


r presents the flat front theory. 





evaporator condenser reserwoir 


perature 





Vapor-gas flat front model 


figure 2 
2 





More recently, Marcus and Edwards [Ref. 7] have developed 
тэг 121056 front theory. Іп this analysis the contribution 
of binary diffusion between the non-condensable gas and the 
working vapor is E E adompaeuwsbh axial conduction 
YY the pipe wall. In addition to heat transfer from the 
uuenser by convection the contribution due to radiation 
is included. The inclusion of these phenomenon results 
inthe vapor-gas interface spreading out Into a more 
realistic profile. This analysis develops accurate 
temperature and concentration profiles for the condenser 
ew Onl. 

This analysis was used to formulate the TRW computer 
program Шин 01-15 Capanle Of сале а та ЙГ ГР varrous gas 
Neg heat pipe geometries, among other parameters: 

(а) the wall temperature profile; 
(b) heat and mass transfer; 


(c) gas reservoir size. 


C. OBJECTIVE 

When a container is filled with two vapors of dissimilar 
vee cular weight the vapor of larger molecular weight will 
st le to the bottom portion of the container due to gravity. 
In the operating gas loaded heat pipe the non-condensable 
gas tends to collect at the end of the heat pipe away from 
mics evaporator, having been convected there alone with the 
working vapor. An inherently unstable situation will 


therefore be present in the gas loaded heat pipe if the 


Шо 





neenser is above the evaporator and the molecular weight 
The non-condensable pas exceeds that of the working 
Or ие tendeney ol the gravitational force is to 
oppose the convection of the non-condensable gas, thereby 
spreading out the concentration of the gas. This decrease 
in the gradient of the working vapor concentration is 
accompanied by a decrease in the temperature gradient. 
Therefore, it was the objective of this work to investigate 
he operation of a gas loaded heat pipe with particular 
emphasis on the effects of gravity. 

The heat pipe used was designed to operate horizontally, 
vertically with the condenser above the evaporator or any 
angle in between these two positions. This variation of 
the angle of operation allowed a variation in the axial 
component. of the gravitational force. The results presented 
herein are the product of both data obtained from a heat 
J oO designed and constructed Гог this purpose and a one 


dimensional analysis. 
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Iu ED (T LON Or EXPERIMENT 


A. OVERALL DESIGN 

Ber mentioned Ee: ERE а № сопс исенаошко? the heat 
pipe wall tends to spread out the gas vapor interface. То 
reduce this, the heat pipe was constructed from a stainless 
Steel pipe with a wall thickness of 0.0635 inches. Fins in 
the form of Acme threads were cut on a lathe along the 
Condenser section to a depth of 0.0255 inches. 

One of the two end plates was soldered in place while 
the other was designed to be easily removed from the condenser 
mE OI the heat pipe. The removable end plage, as shown in 
Figure 3, was in the form of an internally threaded cap. The cap 
Was fastened to a threaded flange which was Welded to the end 
of the heat pipe. An O ring was placed in the flange to 
ENS CL a pressure and vacuum seal. The eonkenser end cap was 
designed with a 0.25 ineh diameter tube soldered to its 
bernal side. Through a T connection the tube provided the 
necessary inlet for a pressure transducer, vacuum valve and 
Means Of evacuating and filling the heat pipe. 

A mesh stainless steel screen around the inside wall of 
the pipe was used for the wick. To retain the wick firmly 
against the wall a longitudinal stainless steel spring was 
inserted along the entire length of the pipe. 

The evaporator heating device was a resistance heater 


made of .125 inch wide, flat Nichrome strip powered by a 


15 





direct current power supply. Figure 4 shows the overall 


A ign of the heat pipe. 


Outside 
thread 





[DET qu 


End cap design 
figure 3 
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B. DESIGN CONSIDERATIONS 

me performance of a heat pipe is affect@ by any 
@memical reaction or decomposition of the wotking fluid or 
mE corrosion or erosion of the wick, spring, or heat pipe. 
An example is the reaction of water and stainless steel 
resulting ín the release of hydrogen gas. im designing a 
heat pipe the materials used for the pipe and the wick must 
be chosen to be compatible with each other amd the working 
Ве. For this reason, the screen mesh, pile and spring 
Ber were all made of the same material, stainless steel. 

Because the heat pipe operation is dependent on capillary 
pumping the working fluids surface tension head to be consid- 
ered. In addition, the gravitational force &ndà flow losses 
made the consideration of the working fluid's density and 
viscosity important. Methanol, whose properties are slightly 
аг i esirable than water was chosen because of its good 
МИБ characteristic and high compatibili with stainless 
звее!. Ammonia was rejected because of its ttoxicity. 

Following Marcus [Ref. 7], the screen messh size and 
Gamensions of the pipe for a given working fluid can be 
related by a pressure force balance equation. This equation 
ШО Sed on the principle that the net capiigary head in the 
wick must be greater than the losses due to pressure drop in 
the vapor and liquid plus any losses due to body forces. 


This force balance equation is expressed as: 
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э “Аа” мэ... мие А 2 


met body liquid vapor 
capillary - force > pressure = pressure 
head head drop drop 


ho 
where, АР, m 
AP, А Zn Ri 
ЛАШ 
шиг 2 2 
201 Or, ice — R. ) 
2 
апа Шалан 
1 


and the vapor side pressure drop is assumed to be much less 
than the other pressure components and therefore not 
eenendered. 

The heat load was anticipated to be in the range of zero 
to 100 watts with the most probable loading between 10 and 50 
watts. The stainless steel screen mesh was available in 
various sizes up to 200 wires per inch. Using this infor- 
mation and the above force balance equation an iterative 
design ensued for an operating temperature of lo N, 
Various sereen wire spacings, lengths of pipe and number e 
layers of mesh were combined with various piresi teres И 
seeking the proper design. The design centered on providing 
for a pipe from one to two feet long and with a diameter in 


excess of one inch if possible. 
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With a margin of safety, the design was tentatively set 
for a two-inch diameter pipe of overall length 15 inches with 
5 layers of 150 mesh screen. This length allowed for a ! 
Ше) heater section and a one inch adiabatic section. This 
was the tentative design because these parameters were for 
a heat pipe not loaded with a non-condensable gas. To 
achieve an optimum design for the gas loaded heat pipe ihe 
design program developed by Marcus [Ref. 7] was employed. 

In order to use this program the non-condensable gas had 
EUM Selected. The choice centered around the inert gases 
since they were chemically inactive and provided a good 
peSecrum of molecular weights. To assist in emphasizing the 
effect of gravity on the vapor-gas mixture a large difference 
Пишите molecular weights of the working fluid and the non- 
condensable gas was needed. The choice of krypton allowed 
a molecular weight ratio of krypton to methanol of approxi- 
mately 2.5. 

With the heat pipe parameters tentatively determined, 
Marcus' computer program was used to first determine the 
additional length of pipe needed for the non-condensable gas 
шиг ben to give a performance prediction. Two inches were 
added to the condenser to allow for the proper performance 
of the heat pipe with the non-condensable gas. A check of 
the preceding force balance equation used in the design 
process indicated the heat pipe would still operate satis- 


factorily with the extra 2 inches of length. 
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Ge INSTRUMENTATION 

The instrumentation of the heat pipe included both 
Temperature and pressure sensors. Seventeen thermocouples 
memes placed in line in the grooves of the fins at an interval 
of approximately 0.650 inches on the top of the pipe. An 
@o@rtvional thermocouple was placed inside the pipe. This 
Миг шосоир те, enclosed in a stainless steel sheath, was 
Higserted 2 inches through the evaporator end plate. The 
condenser thermocouples were 0.010 inches in diameter copper 
and constantan wires welded with a bead at the end. The 
bead was made with a Dynatch thermocouple welder. These 
beads were attached to the grooves of the condenser with a 
Unitek model 1065 welder. The other end of all the thermo- 
Couple wires vere led to a Jones strip and from there to a 
Hewlett Packard 2010C Data Acquisition со, 

Several layers of Johns Manville MinK insulation were 
в peda around the heater to minimize heat losses. Two 
additional thermocouples were then placed across a layer of 
insulation. By knowing the temperature difference across a 
layer of insulation as well as the thermal conductivity and 
thickness of the insulation the heat loss from the evaporator 
could be estimated. - 

To determine the power input to the heater accurately a 
known resistance was placed in series with the Nichrome 
Meter strip. Ву measuring the voltage aerocss the known 
nccacror the heater current was determined: This. current 


was used with the voltage reading across the heater to 
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determine the power input at the evaporator section. A 
Hewlett Packard digital voltmeter was used for these 
measurements. 

The pressure was measured by a Pace KP-15 transducer and 
aean aed On a Hewlett Packard 1/02B recorder. Figure 5 shows 
a schematic diagram of the heat pipe and its associated 
equipment. Included in this diagram and in Figure 6 is the 
apparatus used to fill the heat pipe with the required 
amount of methanol and non-condensable gas. The filling rig 
used was a previously existing device used by J. Woodard 
[Ref. 8] which was modified to include an inlet for the non- 


condensable gas and an additional vacuum valve. 


D. CALIBRATION AND ERROR ANALYSIS 

The. thermocouples were calibrated in a constant tempera- 
ture bath with a platinum wire resistance thermometer and a 
Rosemont commutating bridge, while the thermocouple readings 
were made on a Leeds and Northrop K-3 precision potentiometer. 
Ine thermocouple readings were found to vary from the com- 
mutating bridge temperature readings by no more than +0.11°F. 
The data acquisition system was then checked against the Leeds 
and Northrop potentiometer and found to have an error bound 
of % 0.49Е, The resistor placed in series with the evaporator 
heater was also calibrated on the Rosemont commutating bridge. 
The resistance was read as 2.013 + 0.0002 ohms. The Harrison 
pressure recorder and Pace transducer were calibrated with a 
Wallace and Tiernan gauge. The error for the transducer and 


recorder was + 0.4 psia. 
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An error analysis was performed for all the data 
recorded as well as the other parameters used in computations. 
A summary of the error bounds is shown in Appendix A along 
Sample calculations. The data points on the graphs were 


drawn to include the error bounds of the information presented. 


E. EXPERIMENTAL PROCEDURE 

Prior to proceeding with the experiment, the amount of 
Pomme fluid to be used was determined. The equivalent 
ШЕ volume of the working vapor for a heat input of 
20 14005 was added to the void volume of the wick structure. 
The volume of methanol required vas approximately 45 milli- 
K depending on the safety factor used. Operation of 
ШИЕ pipe with 30, 40, and finally 55 milliliters of 
mol produced no discernible difference in the temperature 
and pressure readings. Since wick dry out, caused by a 
shortage Of working Fluid in the wick эс, cd" cause 
a breakdown in the capillary pumping process the choice was 
nu ouse 55 milliliters of methanol. An excess of 
uus fluid would not produce any detrimental effects on 
the heat pipe's performance. The krypton gas used in the 
EPIment was research grade while the methanol working 
ша was speetrophotometric quality. 

ог uo assembling the heat pipe, the pipe, wick, end 
Cap, and spring were all scrubbed with acetone and rinsed 
three times with distilled water. The assembled heat pipe 
was subjected to the same process prior to attaching the end 


plate. The procedure used in preparing the assembled heat 
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НЕ (ог орегафјоп started with evacuation of the heat pipe 
amet tlling assembly. The filling assembly, in addition 
tozsupplying the working fluid and non-condensable gas, 
provided for the connection of the heat pipe to the evacua- 
F BUMp. A total of approximately 30 inches of 0.25 

inch inside diameter vacuum hose connected the heat pipe 
through the filling assembly to the evacuation pump. The 
prep was able to achieve vacuum, read at the pump, of 

ж 10-6 millimeters of mercury. 

When the evacuation vas completed the heat pipe and 
filling assembly were isolated from the evacuation pump. 
measured amount of methanol was then released into the 
heat pipe. The pressure was recorded vhen the temperature 
mede the pipe equalled ambient temperature . The krypton 
was next released from its high pressure cylinder and allowed 
tO flow into the heat pipe, while monitoring the pressure 
and internal temperature. The mount of kngypton present 
ша then be calculated by using the temperature and the 
difference between the pressure reading and the partial 
pressure of methanol with the perfect gas equation. 

Wen the internal temperature of the pipe was equal to the 
ambient temperature, it was recorded along with the pressure 
in the heat pipe. 

The heat pipe was then isolated from the filling assembly 
by means of the vacuum valve on the stem of the heat pipe's 
end cap. The pipe was then positioned for the desired angle 


of operation and the heat input was applied. 
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ites temMperavure ana pressure readings were monitored 
шта: а steady state operation was detected. The steady 
state condition was based on the condenser temperatures 
rate of change being 127 Chan -0.9 SR per hour. 

The temperatures and pressure were recorded at steady 
шише then the heat pipe was either repositioned or the 
heat input was altered to initiate a new operating mode. 
This process was continued until sufficient data had been 
obtained. 

Ehurdown of the pipe required only the removal of the 
heat input and subsequent release of the internal fluids 


when the heat pipe had achieved ambient temperature. 
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ЕЕ. ЕЕ Е ТА RESULTS 


A. CONVENTIONAL HEAT PIPE 

The heat pipe was initially operated as a conventional 
75e With only methanol present. Heat inputs were 
s ed from 10 to 50 watts in 10 watt increments in each 
E NHPee heat pipe positions; horizontal, cemdenser raised 
45°, and vertical. The data presented in Figures 7, 8, and 
9 is representative of all the condenser temperature readings 
ET The various heat inputs and positions mentioned. As 
these Figures show, the entire condenser operated in an 
isothermal manner as expected. 

Figures 10 and 11 show the evaporator temperature and 
K Ser pressure for a particular пеас input. The data 
amps SO obtained for each of the heat pipe positions are 
connected with smooth curves. The evaporator temperatures 
were used to determine the saturated vapor pressure of 
methanol from the data on methanol presented by J.M. Smith 
[Ref. 9]. The pressures obtained from J. M. Smith's table 
were all within + 5% of the corresponding pressure readings 
One E ге l0. Therefore, the assumption of operation of 
the heat pipe with a saturated working vapor is confirmed. 
The temperature readings used to enter J. M. Smith's table 
were from the evaporator while the pressure reading used 
for comparison were taken from the end of ра сопдепзег. 


The fact that the pressures are equivalent, within the 
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e@emracy of this experiment, and were read at opposite ends 

o61 the pipe indicates that there is no pressure gradient 

throughout the conventional heat pipe during operation. 
Eure Ji indicates the heat input to the heat pipe 

Беше tion oí the evaporator temperature for the herizen- 

tal and vertical positions. The variation of the evaporator 

vemperature for the changes in heat input is shown to be 

linear. The slope for the vertical position in Figure 11 

me, Sess than that for the herizontal position. This indi- 

cates higher evaporator temperatures for the same heat 

input in the vertical position when compared to the horizon- 

ог 00511101, An explanation for this can be found in the 

heat transfer coefficient of the cylinder which is directly 


ва Со the Nusselt Number by: 


where the Nusselt Number equals: 


T25 
horizontal Net, = 525 (МР, Nay? 
225 


vertical Non 59 (NS, Nap? 


where: 
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Lure L is the significant length and equals: 


Hn pipe rameter Tor the horizontal position, 


где оуегайт ріре length Tor the vertical position. 


Since the overall length of the heat pipe is approximately 
mes its diameter the heat transfer coefficient for the 
horizontal pipe exceeds that for the vertical case. For the 
same heat input the decreased heat transfer coefficient 
ue vertical position, when compared to the horizontal 
Eu lon, requires a greater temperature difference between 
the heat pipe and the atmosphere. This is reflected in 


Figure ll. 


B. GAS LOADED HEAT PIPE 

After evacuation the heat pipe was filled again with 
НИ 1 11 Сетз ог methanol and 2.5 x 1073 poundmass of 
Krypton gas. The heat pipe was operated over a similar 
ие of heat inputs and the same three positions. The 
condenser temperature profile for various heat inputs and 
positions are shown in Figures 12, 13 and 14. 

Figures 12 and 13 representing the horizontal and 45° 
position, show only a nomial negative temperature slope 
from the evaporator end of the condenser. This temperature 
profile indicated the concentration of working vapor was 
decreasing very little in the direction of decreasing 


ш шрегасите. it was anticipated that a larger temperature 


gradient would be present, reflecting a more marked 
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segregation of the working vapor and the gas. It was further 
алсаз расеа that this front between the vapor and gas 

wouUutd spread out over more of the condenser as the angle 

Of the heat pipe increased. 

The vertical operating position, as shown in Figure 14 
had a larger negative temperature slope than that found in 
Che ether two positions. The vertical position appears 
to indicate the absence of a sharp front but rather a 
meure of the vapor and gas. The concentration of the 
vapor does decrease more markedly with distance from the 
ао? in this position. This indicates that inspite 
©f the gravitational force and the larger molecular weight 
Ae he non-condensable gas, the gas does tend to accumulate 
in the condenser at the end away from the evaporator. The 
concentration of the non-condensable gas at no point excludes 
puerco ndaensatíion of the working fluid in the condenser. 
Therefore, the gravitational force does effect the variable 
conductance property of the gas loaded heat pipe. 

In comparing Figures 7, 8, and" 9 with Pres 12, 13 
and 14 it is apparent that the change in temperature of the 
condenser for a variation in the heat load was less for the 
gas loaded heat pipe. This is also indicated by the larger 
Slope found in Figure 11 when compared with Figure 16. This 
"mother feature of the gas loaded heat pipe due to its 
weerabie conductance. As these two figures indicate, the 
operating temperature of the evaporator is higher, for the 


same heat input, for the gas loaded heat pipe. Figures 10 
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(о indicate that, as expected, the pressure is greater 
Іш (пе gas loaded heat pipe compared to the conventional 
heat pipe for the same heat input. 

Ene laek of a Sense temperature variation in the 
condenser for the horizontal position vas investigated 
NENNT after the data had been collected fem all heat 
вис and positions. 

The heat pipe, in the horizontal position was rotated 
атома jts longitudinal axis in order to allow the thermo- 
couples to take readings along the bottom of the heat pipe. 
In this position the readings were uniform and on the 
average 12 °F lower than those taken on the top of the 
в моторе. The bottom of the pipe was expected to be cooler 
Ene top due to the natural convection. “However, the 
temperature slope on the top of the pipe and the absence of 
such a slope on the bottom indicated that the non-condensable 
gas and methanol may have stratified. That is, the krypton 
may have settled on the bottom of the entire length of the 
pipe”. 

Шис пе gas vapor had not stratified the nominal tempera- 
ture gradient present on the top of the pipe would have 
been present on the bottom of the pipe. With stratifica- 
tron, the heavier krypton would settle to the bottom portion 
ОТ (пе horizontal condenser. This would greatly reduce 
the condensation rate if not completely stop it in this 


Section. The bottom of the heat pipe being greater than 
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ambient temperature can be accounted for by the circumferen- 
tial conduction of the pipe. 

To further @ivesuarate this, temperature sensitive 
meta crystals werc applica over a 2 inch wide black paint 
o Which had been applica vo the entire length of one 
side of the condenser. The liquid crystal has the charac- 
Цио лес of changing Irom clear to a color over a small 
temperature range. The crystals used were R53. These 
indicated Спргее ари те temperature ranges by three 


пи ont colors as indicated below [Ref. 10]: 


Color Temperature °C 
Red DOS ретро но 

Green Den. TOES 
Blue а ыб 


With the liquid crystals on the pipe the hedik input 
was varied until the general temperature range for crystal 
coloring was achieved. Figure 17 shows two photographs of 
the isotherms observed for the indicated power input. Іп 
the case of a heat pipe the isotherm represents a line of 
constant working vapor density. Dievieceusrmat Fhe isotherm 
lines run down the side of the pipe in a near horizontal 
manner indicates the working vapor density varies in the 
same manner. Stratification to some degree may therefore 


be present. 
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Heat input equals 21.70 watts 





Heat input equals 23.6 watts 


Photographs of liquid crystals on condenser 
figure 17 
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The effect of natural convection upon the slope of the 
iostherms has not been fully Invest 2 02:t chis time and 
would provide more insight into this phenomenon. The 
m nuiseastjions of the possibility of stratification will be 


СЕЕ СИБ сеа 1л Section V. 
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IV. ANALYSIS 


ШІ THEORY 

The analysis developed in this section is for the con- 
denser section of a gas loaded heat pipe. The development 
mS Tough the application of the continuity and momentum 
equations for one dimensional flow, the energy equation for 
EN all, wich, and fin heat transfer as well as Fick's 
law of diffusion and the perfect gas lav. The development 
Бог 0160 111 four first order differential equations which 
Е Solved simultaneously with the use of a Runge Kutta 
poutine. The formulation of this analysis is outlined 
below. It should be noted that as Figure 18 indicates, the 
ЕП Of the coordinate system is at the end of the pipe 


away from the evaporator. 





P^ 


Heat Pipe 
figure 18 
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p ji with the integral form of the mementum equation 


applied to the section shown in Figure 19: 





Р %у 
р _ 
п e 
ХУ шаг 
— 7 ك‎ 
252 — 
E а 
а (ФУ) — 98 
+ ах = 
+ AP 
4 ах 


figure 19 


e» JF, = Í oV(V:dA ) 


In the following analysis the viscous term has been left 
out of the above equation. This simplifying assumption is 
EI gon the Reynolds number being very small in the 
operating range of the heat pipe. For a typical heat load 


ОГ 20 watts: 


AUN a 
Rp с Криста 5020 
Q 
Where, V = 
Dep А, О 
ү = 68126 
B T 2 
(181) 7 (1/6)° (1/20) 
12-12 s Seo 
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| Го ју developed laminar flow the friction factor f 


can be represented by 64/Rp . This results in a shear force: 
Е = буту dx 


што деглпе the entire one foot length of condenser 
the force is in the range of 1079 pounds or less. Therefore, 
the shear force can be neglected. 


Referring to Figure (19), equation (1) becomes: 


ВИМЕ СЕ ах] + ок 2јп 6 dx » -pV^ + [oY + A (oV) ax] 
dx ах 
ay . 
T e ” 
[V + ES dx] + m ps т Dido 


where P is the pressure in the condenser section, V is the 
vapor-gas velocity, Ма is the vemocjty of the Vapor, р 155 
Napor=sas density, g is the gravitational acceleration and 
D is the diameter of the vapor space. 

Dividing by A,dx, and taking the limit as dx approaches 


Zero saves: 


ЧА ар 1 = dV alev) 1 мас 
(2) за + pg sin 0 pV UN V dx + m us 


>|: 


mee cONCinuity equation for the mixture is: 
27: ЗЕ 
рУА, = mvm D dx + [pV + zs CON dx] А. 


OI 


e) По 
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Applying Fick's First Law to species a: 


- NR 
(4) Ба DE 


where 89) Ts the diffusion coefficient for the vapor-gas 
mixture, p, is the density of the vapor and w is P/P. 


An alternate form of this equation is: 


_ АЯ dw 
(5) V(w-1) = DY 


The pressure in the heat pipe condenser, P, is the sum 
of the partial pressure of the vapor, Бар amd the partial 


осе of the gas, Ры: 


The perfect gas law can be used for the non-condensable 
eS. This is possible because the operating temperature and 
pressure in the heat pipe are well below the critical point 
of the gas. However, during heat pipe operation the working 
1110 уарог is in a saturated state which precludes use of 
the perfect gas equation. Instead the pressure of the 
working fluid vapor can be represented as a function of 
temperature only. Therefore, - 


R 


(6) роди LU. M. T 

ог, 

(T) P= f(T) +g- p (1-w)T 
b 


ш 





те В is the universal gas constant, М, is the molecular 

weight of the gas and f(T) is an equation based on the 

saturated vapor data presented by J.M. Smith [Ref. 9]. 
Sensaider the section of wick-wall shown in Figure 20 for 


ME nergy balance. 


| TACK 
i ЇР 
иа аи с. 
h (T- T4.) all 
fme 220 


Assuming the temperature difference between the wick on 
the inside of the pipe and the outside of the pipe is negli- 


Беси the radial direction: 


ПЕ m he TD dx = q + dq + h(P-T,) пор dx 


ШЕ 


here. q is the heat being conducted axially by the pipe, 
m is the rate of condensation of the woking рег иас 


area, h is the latent heat of condensation of the working. 


fg 
mor, h is the heat transfer coefficient, T is The ambient 


temperature and Do is the outside diameter of the heat pipe. 
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lom Fourier's conduction law: 


(9) 20202 - 


where E is the eross seetional area of the wick-wall and 
Нэгэ їе thermal conductivity in the axial direction of the 
"well combination. The evaluation of the thermal 
@emevetivity,-k, and the heat transfer coefficient are 
шешеторей ас а later point in this section. 

Substituting eauation (9) into equation (8): 
б, ат. h ER A 4 m Dep e D 


S 
ае wo’ L m k Ay 


The equations necessary for the development of the 


analysis are summarized below: 








ар ë - оу 4У а(оу) ртр у 
(2) ео ах A. кт 
° TD d 
(3) m NW - ~ ax (oV) 
C 
E dw 
(5) у= у ша 
(6) P = f(T) + = (l-w)? 
b 
2 h тр 
EI eT... e {ү 5 
we k A ГА 
X p p 
апа, 
(11) p = 5 


EL 





DEUS ia an equation for p, as a зао о Cie! 2 based оп 
the data presented by J.M: Smith [Ref. 91. 

СОСЫ Оеп тт une above equations, the diffusion coerfi- 
E D , has been taken as constant. The diffusion 
ШОО C ent equation for low pressure is based on the 


wine equation [Ref. 11): 


1.858 10-3 11.2 [ pe ] 
9 рс: в 
2 
о р 


lungs indicates that) does vary somewhat with temperature. 
eT, the above assumption has been made in order to 
ту the analysis. A mean diffusion coefficient for the 
шпрегасиге гапре involved was used in the program. 

The following non-dimensional variables were used in the 


above set of equations: 


p* = 5- 
O 

x 
о" = L 
O 
x Sx 
ШЕГІ 

х 
122 г- 
р 
m, 
Мете” ше = — 
Ро 

: Q 
and, m = 
5 Г Dep Ас 
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m 
m = — h+ƏÜ> 
m 
г 
E X 
x = f 
E D 
аг р 





‚ Р х 
ТЭГ... | 2122 r — 11.) 
dx 





В пре р Nee gb 
k sx x 
dx dx 
HQ 2: Мез T құ 
TD p hog Š 
Es х + 
(13) т = - үт в (ру) 
EET a2 " D^ p, пр o, he 2 au 


SD ATL © m 


x R Po Ша х ж 
EORR e D -1--2-21 T (I w) 
M Po 
* 3 
аёт E hL A. * 1 Ц 1,2 Q es 
ns TA И 
dx p о Dk A, T, 


substituting for the non-dimensional coefficients in the 


preceding equations: 
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gL 
T D ZI 
ЦТ, 
ЕС 
"UR 
Ц КА 
р 
4 co 
5 Dk А Да 
3 c 
д пр робе AD) 
6 169L° 
N pita Polo 
7 M, Ва 
прер в, 20) 
N.N Ou. 
227 QL 
the equations become: 
: 5 
(17) р sin 6 - N, => = М> [o V dV. o. ve dio ) 1 
ых X ах 
. x x 
(18) m = - Ф (pV) 
Зах 
(19) ү* ES 1 N.N ди 
е SNO dx 
к x х x 
(20) ЈЕ = EUMD EE Ny (1 + w)P T 
aer” = + OK 
(229) = № T T )- Nm 
dx 2 1 D 2 
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Investigating the relative magnitude of the coefficient 
of the momentum term in equation (17) with respect to the 
E ceoefficients for representative values of the operating 


moat pipe, one finds: 


ТЕ 4x10x3. 413 Jr 1 1 
2 ШОСТЕ ОТИ 1 25348081 (32.2)(3600) 
Е =, 
N, = ЦО 


À similar investigation of the other terms indicated 
that Nə was at least a factor of 107? times smaller than the 
Нэг 101 Lhe other non-dimensional coefficients. This led to 
EN pliifylng assumption that the momentum term could be 
excluded from the analysis without loss of accuracy. The 
emulation of this coefficient Was retained in the final 
program in order to serve as a check on its relative value. 


The momentum equation then reduces to a force balance: 











x 
(22) p sin @ - N, “y= 0 
ax 
Differentiating equation (20): 
ар“ аг (т) aT. И < as 
Е NAS + пана 
dx ат dX ат" dub 
1 «Ка? dw 
= зи pr POE 
T w* ах 





Subtracting this from equation (22072 


x % * 
e * Y 
(24) = Š sin Ө T (1 (8 + 99, те) = 
1 dT ах ат ах 
хх 
E Ny = а m Aue, 
и dx 


Solving equation (2H) for EN and differentiating: 


ах 
2 2m“ 
72159 A 136. + w(l-w)G,] om 
2 17952 Ц 2 
ах 7 ах 
2 ат 2 
+ < G. + vw(1- w)65JC—g) 
7 dx 
2 ат’ а 
+ [° Z G. + (1 - 20w)G, J = = 
N. "1 Li : 
7 ах di 
+ Po pz l sin 22-21 —5 sin 0 a] 
177 гр ах T ах 
nere, 
х 
к 
| ~ пе КЕ 
desi. p 
2 3 5 
іе се (05205205 dM. 
ү - cu cw - : 
24:01 аат ат 
G4 = z 
Coa) š 
1 as” 2 
бу, ЕЕ F 
d (118 о 
x x 
ЕЕ 45 2 
жо NEED! 
ТЫС aT таз ат ат 
5 S e 
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Swbetitutinme@eequation (18) into equation (21): 


2-1 N 
d х X = х х 
(26) Пор = (Т - т, ) + ро Y) 
ах | ДУ 


Differentiating equation (19) and substituting into 


equation (26): 


2 х х | x 
eT х х 95 ji aw dT 
T. 2 lj 5 ar’ ^ w-1 a 

гэ 3 је (67 — 1) E 


ENSStituting equation (25) into equatd3n (27): 








M . f x 
dus. * * dT dw 
(28) Gio С ято. = U - 1) (7 = а ) + 044 — е 
dx dx dx 
170 aT". 
UM o 
le ах 13 ах 
E 
- Gay Sin 9 - + Gis Sin 6 сш. 
dx dx 
where, 
N-N 
Р 56 % x 
Gio Kr Ny wS G) + (w - 1)N¿N¿S Gy J 
x 
G = 56 dS TE $ с + N_N E С 1 
11 W Е N 1 5 6 W lj 
ar T 
RS x 
ба = МН A 1 с 
w (w-1) 
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NEN 


ты с” | _ Ч 
б; Г Ny wS G + NN; GUESS Gc] 
x 
a, = Oi Ss 
14 NN, ү! T^ 
15 NN. mre 


Sumplifying equation (25): 


x ¿ x 
EE uw... 91 o ат, 92, 9v 
_ #2 20 EROS 0] ww E 22257 8 
ак ах ах ах (Ж 


P 
dw | ат 
- 6 sin ONE x 5 Goy en 00 F 


23 dx dx 
mere, 
we 
Gap ® E E w(1-w)8,] 
í 
ме 5 
854 = If G3 + w(1-w)Gg] 
7 
ом 
G55 - Ын G; + (22-27 10) 
Ї iJ 
G = [—— —]J 
23 № т 
1 W 
رر‎ = I 5) 
21 N. Nç Т 2 
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Bemeations (20) and (29), along with their associated 
Memeaamensional coefficients, form a set of two, second 
Order, simultaneous non-linear differential equations. As 
mentioned earlier these equations were solved using the 
IBM 360/67 computer with a Runge Kutta routine. The entire 
program was put in double precision to avoid as much machine 
ONA off error as possible. 

NI value of the heat transfer coefficieme for natural 
convection, h, was calculated based on the Nusselt Number, 
a лег. 12]. 

GE 


where, 5 mn cp ов 


where, m .59 for the "eetical 510101 


and, m ¿Bes for the horizontal position 

The Nusselt Number for the heat pipe in the horizontal 
pon was found to be 11.59, while for tWe vertical 
position the Nusselt Number equalled 67.84. These values 
Nererused in the following equation to fina the natural 
eae ection portion of the heat transfer coefficient, Пс: 


To this was added the value of the radiation component 
Of the heat transfer coefficient, Ше restmreaneydn the 


following final equation: 
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OE о 


D DM р =) 


Two values for a condenser temperature of 600? R and an 
ambient temperature of 540°R are indicative of the relative 


values used: 


horizontal position h = 1.41 


vertical poSI ion h = 1.00 


de terminine the axial thermal conductivity of the 
Ш ehd saturated wick combination, k, the conductivity of 
pou asvurated wick was first determined and then combined 
КИЕ И The thermal conductivity of the stainless steel wall 
Miwa perallel manner. The determination of the saturated 
wick thermal conductivity could be viewed as either a 
parallel or series combination of the wire screen and the 


working fluid [Ref. 7]. The two values obtained are: 


| ФЕ 
k - panl = „296 
series ФК. + (1= 6) Kr у 
Koarallel y (1-9)k, a = 5.80 


Mere 9 15 the*voWd fraetion in the*wYiek, and is eqwual 
во. 397. 

In assuming the value is that of the parallel case and 
combining this with the thermal conductivity of the stainless 
steel wall, the largest possible value of thermal conductivity 


Tor сле WICK wall Combination is obtained. This is what the 
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no 


Exe of the fins was intended to reduce. The axial 
conductivity used was therefore: 


= ВТУ 
о es 


БОО RESULTS 

Computer program developed to solve the set ог 
simultaneous, first order, non-linear differential equations 
is included after Appendix A. Several solutions obtained 
from this program are included in Figures 21 and 22. The 
Solutions vere obtained for a specific heat input by assuming 
Р erature for the condenser at x equals zero. As a 
first estimate the condenser temperature from the experimen- 
шештетеве for the same heat input was used. The computer 
s Op indicates in the output the amount of heat the pipe 
"ut pansfer at this temperature divided by the amount of 
heat specified as being the input. The condenser temperature 
Ага until this quotient equalled one. 

Meer analytic solutions shown in Figures 21 and 22 
indicate a nearly isothermal condenser. The condenser- 
evaporator interface is only from 0.01 $0 0.2% greater in 
temperature than the other end of the condenser in these 


figures. 
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V. DISCUSSION 


A. COMPARISON OF ANALYSIS AND EXPERIMENT 

Figures 12 and 14 present condenser temperature profiles 
based on the operation of the gas loaded heat pipe. Figures 
21 and 22 present the condenser temperature profiles based 
on the analysis developed in the preceding section. The 
B urrIons and heat inputs are the same for Figures 12 and 21 
and the same for Figures 14 and 22. 

The analytic solutions all have higher condenser 
temperatures than the data obtained from the operating heat 
mire or identical heat input and position. The analytic 
condenser temperatures exceed the operational temperature 
in °F by 7.2 to 19 percent. The larger difference between 
Miemanalytic and operational condenser temperature is 
вона 2n the vertical operating position. 

Шие heat transfer coefficient, h, if increased would 
educe the condenser temperature since for a given heat 


ораистег, Q: 


Q-n(T-T7) 


The values of h determined earlier were based on natural 
Son vection over a symmetric cylinder. The addition of fins 
and the removable end plate with its associated flange on 


the pipe served to disturb the flow over the pipe. This 
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Пеш оапсе ог the natural convection over the heat pipe 
would have more effect in the vertical position where the 
s tion flow travels the length of the pipe. Any such 
disturbance would tend to allow greater mixing in the con- 
vective flow. This would allow a greater amount of heat 

to be transferred at the same condenser temperature or the 
Original amount of heat transferred at a slightly reduced 
(ита сиге. This is the effect of increasing the value 

її the preceding equation. Increasing the value of 
the heat transfer coefficient used in the analytic solution 
will reduce the condenser temperature obtained. This 
DEAS (о therefore be one of the reasons for the differ- 
e in the analytic and operational condenser temperatures. 

ME analytic condenser temperature was nearly isothermal 
эш earlier. This again differs from the operational 
condenser temperatures obtained. The condenser temperature 
ES Ses vith increasing distance from the evaporator іп 
wm The analytic and operational solutions. However, the 
ШОШО ЛО solution did not have as large a decrease as the 
operating heat pipe did for the vertical position. 

The difference in the temperature slopes of the analytic 
solution and operating heat pipe are noe of such magnitude 
momtrequire the analytic solution to be discarded. In fact, 
Boch solutions indicate the absence of any front. Therefore, 
the assumptions of negligible viscous and momentum forces 


in the analysis are justified. 
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Ee analysis is unable to give any insight into the 
рама пра 11су of stratification in a horizontal gas loaded 
Meat. pipe, except to confirm the lack of any vapor gas 
meno. tne possibility of stratification in an operating 
heat pipe cannot be overlooked. The design of' the gas 
reservoir may require alteration to allow the stratified 
ИСО flow in and out with heat input and pressure 
changes. 

The heat transfer characteristic of a gas loaded heat 
uu bh stratification will differ from the flat front 
ствие models discussed earlier. The stratified heat 
pipe may produce some of the characteristics of the 
diffuse front model with the front turned 90° in the 


pipe. 


B. RECOMMENDATIONS 

Further investigation with the computer program has 
шог ted a heat pipe of different dimensions from that 
used in this experiment may yield more interesting results. 
A IOonger pipe of smaller diameter with a higher heat trans- 
fer coefficient produces a marked vapor-gas front in the 
EM tro solution. With such a diffuse front the rotation 
of the pipe through 90° could be 2 to observe the 
gravitational effect on the vapor-gas front. 

addition, other inert gases could be introduced into 
рас ехіссіпс heat pipe in place of krypton. Lighter gases 


may allow the vapor-gas front to become more apparent while 
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Bas phenomenon. 


67 





APPENDIX A 


SAMPLE CALCULATIONS OF ERROR ANALYSIS 


2 172 
2 N 2 -V 2 
E Ou ша: 


mere Y is the voltage, R is the resistance, Wp is the 
ГОСТ bound in the resistance and Wy 156 the error bound 
in the voltage. 


fn T (ing both sides by Q: 





W W 172 
Q 252 R. 2 
ioe UG WU) t (q? ] 
W 172 
Q 2 2 003,2 
10,815 5.02) a | 


W 
ES 

б 0.0093 
W 

277 | 
св 0.93 7 
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Quantity 

T Temperature 

Р Peessure 

Q Heat input 

L Length 

D Diameter 

%) Diffusion coefficient 

h Meat transfer coefficient 
E Gravity 

Deo Latent heat of vaporization 
g васе tension 

р Density 

V Kinematic viscosity 

V Voltage 

R Besıstance 

Ши Void fraction of wick 

m Condensation rate 
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